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Abstract 

We consider effects of the fields of strong electromagnetic waves on various characteris- 
tics of quantum processes. After a qualitative discussion of the effects of external fields 
on the energy spectra and angular distributions of the final-state particles as well as on 
the total probabilities of the processes (such as decay rates and total cross sections), 
we present a simple method of calculating the total probabilities of processes with 
production of non-relativistic charged particles. Using nuclear /3-decay as an example, 
we study the weak and strong field limits, as well as the field-induced /3-decay of nuclei 
stable in the absence of the external fields, both in the tunneling and multi-photon 
regimes. We also consider the possibility of accelerating forbidden nuclear /3-decays 
by lifting the forbiddeness due to the interaction of the parent or daughter nuclei with 
the field of a strong electromagnetic wave. It is shown that for currently attainable 
electromagnetic fields all effects on total /3-decay rates are unobservably small. 
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1 Introduction 



Study of quantum processes in intense electromagnetic fields is a very interesting subject. 
Strong external fields can help us to learn more about the properties of the involved par- 
ticles and their interactions. Studying processes in strong fields may also have interesting 
implications for astrophysics and cosmology. Recently, there has been a renewed interest in 
this topic in connection with development of new powerful laser sources. 

In this paper I will discuss effects of strong external electromagnetic fields on various 
characteristics of quantum processes. In sec. 2 a rather general qualitative analysis of these 
effects is given, whereas sees. 3 and 4 are dedicated to a specific example - nuclear /3-decay 
in the field of a strong electromagnetic wave. 

My interest in this topic was raised in the early 1980s by A.M. Dykhne, who called 
my attention to a paper published in the Physical Review Letters [1]. It was claimed 
in that paper that under the influence of electromagnetic fields of existing at that time 
powerful lasers /3-decay of tritium can be significantly accelerated. Simple estimates I made 
did not confirm this conclusion, but at the same time I could not pinpoint a mistake in 
the calculation done in [1]. The problem was that the calculation was very complicated 
and difficult to follow. It was based on the standard at that time procedure of infinite 
summation of partial probabilities corresponding to absorption from the external wave (or 
emission into it) of all possible numbers of photons. This motivated me to look for a simpler 
way of calculation of the total probabilities of quantum processes in the fields of intense 
electromagnetic waves. My quest was strongly supported by Victor Khodel, a colleague of 
mine at the Kurchatov Institute, who used to say that "if there is a simple result, there 
must exist a simple way of obtaining it" . 

Eventually, I found a very simple way of calculating the total probabilities of quantum 
processes with emission of non-relativistic charged particles. The method in particular 
allowed one to easily study all interesting limiting cases - the weak and strong field limits 
as well as the case of the field-induced decay of a particle (or a nucleus) that is stable in 
the absence of the external field, both in the tunneling and multi-photon limits. The results 
were published in [2, 3]. 

While I was working on this subject, several papers appeared [4, 5, 6], where the problem 
of /3-decay in strong electromagnetic fields was re-investigated and it was shown that the 
results of [1] were erroneous. The analysis in [4, 6] was based on the same old summation 
technique, whereas the approach in Voloshin's paper [ ] was close in spirit (though not 
identical) to the one I was developing. In sees. 3 and 4 I will discuss the method of [2, 3] 
using nuclear /3-decay as an example, but it can actually be applied to a much wider class 
of problems. 

Gratifyingly, in all considered limiting cases the results of direct calculations of the 
probability of /3-decay in the field of a strong electromagnetic wave [2, 3] agreed perfectly 
well with my previously made estimates, often even including the numerical coefficients. 
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1 learned a great deal about how to analyze physics problems qualitatively and how to 
make simple estimates from A.B. Migdal, both from my personal interactions with him and 
from his splendid book [7]. It is therefore a great pleasure and honour for me to write this 
article for the special issue of Yadernaya Fizika dedicated to the centennial anniversary of 
A.B. Migdal's birthday. 

2 Qualitative considerations 

Consider quantum processes such as 

o A + B ->■ A + B (scattering) 
o A + B ->■ C + D + . . . (reactions) 
o A ->■ B + C + . . . (decay) 

How can external electromagnetic fields influence these processes? They can 

• Modify the differential characteristics of the process (i.e. the energy spectra and an- 
gular distributions of final-state particles). 

• Affect the total probabilities of the processes, such as total cross sections and decay 
rates. 

• Finally, new channels of reactions or decays, which were not available in the absence 
of the external fields, may open up. 

Let us discuss qualitatively all these possibilities in turn. 

We will now make several assumptions that will be used throughout this paper. First, it 
turns out that the smaller the characteristic energies of the charged particles, the stronger 
the effects of external electromagnetic fields on the processes in which they are involved. 
For this reason we shall consider processes with non-relativistic charged particles. We will 
assume that the system is subjected to the field of a monochromatic electromagnetic wave 
of frequency u and electric field strength E, and that the field is a low-frequency one: 

HuJ^Ef), (1) 

where Eq it the characteristic energy of the process. The low-frequency limit very often also 
means that the wavelength of the photons of the external field c/w is large compared to 
the characteristic length of the process l x (such as, e.g., the atomic size in photo- ionization 
processes or the nuclear radius in nuclear /3-decay): ul x /c <C 1. The assumption that all the 
participating charged particles are non-relativistic allows us to concentrate only on effects 
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of the electric component of the external field, whereas the condition wl x /c <C 1 implies that 
we can adopt the dipole approximation, in which the external electric field can be considered 
as a time-dependent uniform field E(t). Such a field can be conveniently described either 
in the Coulomb gauge 

A"(f,x) = (0,A(t)), E W = 4^T (2) 

with the coordinate-independent vector-potential A(i), or in the scalar gauge 

A"(t,x) = (0(*,x),O), 0(t,x) = -E(t)x, E(t) = -V0. (3) 

In different situations different gauges turn out to be most convenient for calculations; we 
will use the Coulomb gauge in sec. 3 and the scalar gauge in sec. 4 . 



2.1 Differential characteristics 

Strong external fields can modify energy spectra and angular distributions of particles pro- 
duced in scattering, reaction or decay process. This happens due to the absorption by a 
charged particle of photons of the external field (or stimulated emission of photons into it). 

Free on-shell particles cannot absorb or emit photons due to energy-momentum conser- 
vation. However, a particle that undergoes a scattering which changes its momentum, or 
is produced in some process (such as an electron production in /3-decay or emission of an 
electron from an atom due to photo-ionization or electron-impact ionization) can exchange 
photons with the external field. 

Let as a result of some process a non-relativistic particle of charge e and mass m be 
produced, and let its kinetic energy be e < Eq, where Eq it the energy release in the process, 
i.e. the maximum kinetic energy available to the particle under consideration. The particle 
can receive some energy from the external field or give to the field a fraction of its energy. Let 
us estimate the corresponding energy As. The exchange of the energy between the particle 
and the field effectively takes place during a characteristic time of order of the period of 
oscillations of the external field: t c har ~ 1/w (the contribution of an integer number of full 
field periods T = 2h/uj averages to zero). Therefore the momentum that the particle can 
obtain from the field is of order 

bE 

Ak = eE t char , (4) 

where Eq is the amplitude of the electric field strength. For particles with a characteristic 
energy Eq (i- e - with the characteristic momentum k = \/2m£o) we have 

Ak _ eE ° =C , (5) 
k i/2me u 
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This parameter characterizes, in particular, the modification of the angular distribution of 
the produced charged particle. If £ is not too large, during the characteristic time t c h ar the 
particle moves over the distance I ~ vot c h ar = ^/^o/ m t c har- The energy obtained by the 
particle from the external field is just the work of the field on the particle over the distance 
/, which gives 

Ae eE l eE / [2e^ 1 



2£. (6) 

Eq Bq 6q V V TJl 00 y 

This result is only valid assuming that £ < 1; for £ ^> 1 one has to take into account that 
the velocity of the particle increases with time and is no longer equal to its original velocity 
Vq. In this case 

I ~ , W + -%, (7) 
m 2 

where eE /m is the particle's acceleration in the external field. This yields 



As ^ em = ejo, 2eol + eEo l\ ^ (g) 

e e e \\ m u m 2u 2 J 

Alternatively, this result can be obtained from Ae ~ [(ko + AA;) 2 — fcg]/2m and eq. (5). 

Thus, the modification of the differential characteristics of the process (energy spectra 
and angular distributions of the produced particles) is governed by the parameter £ defined 
in eq. (5). For £ > 1 the external field sizeably affects these quantities. The total number 
of photons absorbed from the field (or emitted into the field) for £ < 1 can be estimated as 

In the considered low-frequency limit (1) it can be very large even for not too strong fields, 
when £ is relatively small. 

The above estimates should be taken with some caution, though. The parameter £ 
diverges in the constant-field limit u — > 0; does this mean that Ak and Ae will diverge as 
well? In fact, the above estimates of these quantities were made for the field of a plane 
electromagnetic wave, which in the limit u — > goes into crossed uniform electric and 
magnetic fields of infinite space-time extension. For non-relativistic particles only the electric 
field matters, and in an electric field of infinitely large spatial size the quantity Ae can indeed 
formally become arbitrarily large. One should remember, however, that in reality all fields 
are limited in space and time and, in addition, the distance between the source and detector 
is finite. This leads to a natural cutoff in the expressions for Ak/ko and Ae/e^ in the limit 
uj — > oo: the characteristic time t c har ~ should then be replaced by the smaller between 
the time scale of the field and the particle's time of flight between the source an the detector. 
The characteristic length I in eq. (7) has to be modified accordingly. 

Note that the parameter £ does not contain h, i.e. is a purely classical quantity; therefore, 
the distortion of energy spectra and angular distributions of the final-state particles in 
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external fields is a classical effect even when No ~ 1. In many interesting cases, however, 
and in particular for /3-decay in strong laser fields, one has hu> <C £ 0) so that for not too 
small £ one has iVo ^> 1. In those cases the exchange of energy between the system and 
the external field has a multi-photon nature. At the same time, the number of photons Ni 
absorbed from the external field (or emitted into it) during the process of the formation of 
the emitted charged particle may be small. Indeed, the formation process is characterized 
by the time scale to ~ fc/ £ o] therefore the number of photons absorbed or emitted in the 
course of the particle production process is 

N, -- eE oV 2£ o/ mt o _ o eE ° = 2£. (10) 

This means that N\ ~ (hw/e )N <C N , i.e. that the number of photons absorbed or emitted 
in the course of the production process is small compared to the total number of absorbed 
or emitted photons. The same is true for the energy change: the energy Se obtained by the 
particle from the wave (or given to the wave) during the process of its formation is small 
compared to the total change Ae of the particle's energy: 

5e ~ eE ^2e /mt ~ — Ae < Ae . (11) 

This means that for ftw/eo <C 1 the distortions of the angular distributions and energy 
spectra of charged particles in external electromagnetic fields are mostly formed after their 
production process is already over. 

Thus, we have the following picture of how a quantum process in a low-frequency external 
electromagnetic field occurs. The whole process can be divided into two stages. In the first 
stage, charged particles are produced in a reaction, decay or scattering process; 1 during 
the second stage, the produced particles exchange some energy with the external fields, and 
the observable energy spectra and angular distributions are formed. The two stages are to 
a large extent independent of each other, though the second stage is, of course, impossible 
without the first one. 

From the above picture it follows that under certain conditions the differential cross 
sections of the processes should have a form of a product of two factors: the cross section 
of the process in the absence of the external field and the field-dependent factor describing 
the exchange of energy between the system and the external field. The second factor is 
practically independent of the character of the first stage of the process, i.e. is universal. 

The condition faj <C eq which implies that the production or scattering of particles 
proceeds during times that are much shorter than the oscillation period of the external 
field, allows one to consider the first stage as a sudden perturbation, or "jarring" of the 
system in the presence of the external field. Such a concept was developed and analyzed 

1 Note that a scattering process A + B — > A + B with some momentum transfer can be considered as 
destruction of particles in the initial state and their production in the final state. 
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in [8], where it was pointed out that for many different processes (such as e.g. stimulated 
bremsstrahlung, photo-ionization, inverse Compton scattering , etc.) in the field of an 
external electromagnetic wave the differential cross sections take the form 

( dVt ) n \ d£t / n \ hmoj 2 ) ' ^ ^ 

where (da/dQ)o is the field-free cross section, n is the number of photons exchanged with 
the external field (n > corresponds to absorption of photons from the field and n < to 
their stimulated emission), and J n (z) is the Bessel function. From the well known relation 
Yl^L-oo Jn( z ) = 1 h then follows that the total cross section (i.e. the sum of the partial cross 
sections corresponding to all possible numbers of emitted or absorbed photons) coincides 
with the field-free one. In other words, in this approximation the total probability of the 
process is not modified by the external field, 

In fact, the condition fiw <C Eq which ensures Aq <C iVo is a necessary but not sufficient 
condition for the field- independence of the first factor on the right hand side of eq. (12); for 
this, one would also have to require that the energy exchange between the system and the 
field in the process of the charged particle production, 5e = Nifko, be small compared to 
the characteristic energy of the process e . We will discuss this condition in detail in the 
next subsection. 



2.2 Field effects on total probabilities of the processes 

How can an external electromagnetic field affect the total probabilities (decay rates and 
cross sections) of quantum processes? There are essentially three possibilities: 

(i) the fundamental interaction responsible for the process gets modified, leading to a 
modification of the transition operator; 

(ii) the matrix element of the process is altered due to a modification of the wave functions 
of the involved particles in the external field; 

(ii) the phase-space volume of the process gets changed. 

As an example of the first possibility, consider charged- current weak interaction processes 
(such as e.g. nuclear /3-decay) in an external electromagnetic field. Virtual W bosons which 
mediate these processes can interact with the field. However, each act of photon exchange 
between the W boson and the field would lead to the appearance of an extra iy-boson 
propagator in the amplitude of the process, and therefore would strongly suppress it due 
to the very large mass of the W boson. A possible exception is the case when the field 
frequency is extremely high: hu > my/c 2 . However, the latter possibility would correspond 
to a process with the participation of a very hard 7-quantum, i.e. this would be a completely 
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different process. Yet another possibility is when the external field, though a low frequency 
one, is very strong. It is easy to see, however, that in order for the field to produce a 
noticeable effect, the field strength has to be E > m'^ v c 3 /eh ~ 3.2 x 10 26 V/cm, which is an 
extremely large value. It is not actually clear if such strong fields may exist in nature. 

The possibility (ii) can be realized, e.g., in the case of forbidden nuclear /3-decay, where 
the interaction of nuclei with external electromagnetic fields may change the angular mo- 
mentum of the wave functions of the initial and/or final nuclear states, thus lifting the 
forbiddeness of the /3-transition. This possibility will be discussed in sec. 4. 

Let us now concentrate on the possibility (iii), i.e. on modification of the phase-space 
volume of the process. 



2.2.1 Phase space change 

Assume that a charged particle produced in a quantum process obtains some energy from 
the external field. Can this actually increase the phase-space volume of the process and thus 
influence its total probability? At first sight, this seems to be impossible: Indeed, before the 
particle is produced, the field cannot affect it; after it has been produced, any change of its 
energy cannot affect the production probability. You can put the particle into a capacitor 
or accelerator and accelerate it to a very high speed - this would not modify the production 
probability because the production process is already over. 

However, the above argument is purely classical, as it implies that the production process 
is instantaneous. Quantum mechanics tells us that in reality the production of a particle 
with an energy ~ Eq takes a finite time to ~ h/eQ. This is actually a formation time of the 
particle. The produced charged particle is not pointlike - it is characterized by its de Broglie 
wavelength: 

h H , , 

*b ~ _ = - 7 _. (13) 

The formation time t can be estimated as the time it takes for the particle's de Broglie 
wave to emerge from the source (i.e. the time interval over which the particle moves over a 
distance of order of its de Broglie wavelength): 



, . h I m h 

*° ~ w "° = = w,- (14) 

The energy obtained from the field during the particle's formation time can indeed 
increase the phase space volume of the process and affect its total probability (or cross 
section). The energy obtained after that has no effect on the total probability - it can only 
modify the energy spectra and angular distributions of the emitted particles. As discussed 
in sec. 2.1, the latter effect is purely classical. At the same time, as we have just shown, 
modification of the total probabilities of processes in external electromagnetic fields is an 
inherently quantum effect. 
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Let us now estimate the energy obtained from the field by a particle during its formation 
time. It is given by the work done by the field on the particle over the characteristic distance 
l x equal to the particle's formation length, which should be of order Xp. However, after the 
whole particle's de Broglie wave has emerged from the source, the production process is 
already over, so for our estimate we take the characteristic length l x to be 

I. - Y ■ (15) 
The energy gain Sed that affects the phase space volume of the process is then 

5e D — eE l x = e ^ ^ (16) 

and the modification of the total probability is governed by the parameter 

Se D eE h 
X = = ir • (17) 

Note that this parameter contains h, as expected. The modification of total probabilities of 
quantum processes in the external fields would be substantial provided that x ^ 1- 

Let us now give a slightly different argument for this. Assume that a charged particle 
is produced virtually with the energy e + 5e instead of the energy e dictated by energy 
conservation. Such particle can only exist during a finite time interval To ~ h/Se; after that, 
it must be re-absorbed by its source. However, if during this time interval it receives the 
missing energy 5e from the field, i.e. 



m de 

it gets "license to live", i.e. can become real. To affect the production probability signifi- 
cantly, 5e must be > Eo, which gives 

Se/e ~ X > 1, (19) 

i.e. the same condition as we found before. 

Let us stress that the parameter x is independent of u and therefore does not diverge in 
the constant-field limit uj — > 0, unlike the parameter £ defined in (5). This is an important 
point. In ref. [1] the tritium /3-decay process 

3 H -> 3 He + e" + P e (20) 

was considered (eq = (M3h — M3h c — m e )c 2 = 18.6 keV), and it was claimed that in the field 
of a plane electromagnetic wave the decay rate W is 

W^W (E )-[l + c o e}, (21) 
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where Wo(eq) is the field- free decay rate and Co is a number of order unity. Were this result 
correct, it would mean that one could strongly influence tritium /3-decay even by arbitrarily 
weak fields provided that their frequency is sufficiently small, which is obviously wrong. 
Moreover, in the limit u — > the result in eq. (21) clearly violates unitarity. 

We have pointed out in sec. 2.1 that the parameter £ governs the modification of the 
angular distributions and energy spectra of the particles in the external field, and that the 
divergence of £ in the limit uj — > does not pose any problems: the fact that all fields have 
finite space-time extensions introduces a natural cutoff for very small u. However, one can 
imagine a situation in which the field occupies a very large space-time region, and in general 
there is nothing wrong with the fact that, e.g., the energy gain of a charged particle in a 
constant electric field can be very large if the particle propagates very long distance in the 
field. For the total probability of the process, this argument would not work: the probability 
is always limited by unitarity, and therefore field-induced corrections to it cannot contain 
the field frequency u in the denominator. From our analysis it follows that in relatively 
weak or moderate fields the tritium /3-decay rate should be given by an expression similar 
to (21), but with £ 2 replaced by x 2 ( we will discuss this point in more detail in sec. 2.3). 
Comparing eqs. (5) and (17) we find 

X = (hw/2e )^ (22) 

i.e. for hu <C Eo one has x ^ £• 

Recall that the total number of emitted or absorbed photons Nq and the number N± of 
photons exchanged with the field during the formation of the emitted charged particle are 
related to the parameter £ as No ~ ^(2e /hw) ^> £ and N ~ £. At the same time, the 
corresponding total energy Ae obtained from the field or given to it and the energy 5e o 
exchanged with the field in the course of the particle's production can be expressed as 

Ae ~ iV hu ~ 2£e , 5e D ~ Nhu ~ x^o ■ (23) 

The above estimates of 5e o apply actually to the case of relatively weak external fields, 
when fez? < Eo, i.e. x J$ 1- Let us now consider 

2.2.2 The strong field limit (\ > 1) 

In deriving eqs. (13)-(16) we were assuming that the de Broglie wavelength of the emit- 
ted charged particle is actually fixed by the energy release in the process Eo, i.e. is field- 
independent. In the case of very strong external fields one has to take into account that 
the energy gain in the course of the particle formation can be large compared to Eo, i.e. 
the de Broglie wavelength of the produced particle is in general field dependent. Imme- 
diately after the particle production its distance from the source / is small compared to 
its de Broglie wavelength. As the particle moves away and its distance from the source 
increases, its de Broglie wavelength decreases because the particle's momentum k increases 
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Figure 1: Schematic representation of emission of the particle's de Broglie wave from the 
source. Left: weak field case, the de Broglie wavelength is constant. Right: strong field 
case, the particle gains significant energy and its the de Broglie wavelength decreases as the 
particle moves away from the source. 



as k ~ ho + eEot (see fig. 1 where the particle is depicted by a small blob). The total 
probability of the process is only influenced by the energy that the particle gains over the 
distance of order of its de Broglie wavelength, i.e. when the increasing / and decreasing ^ 
"meet" each other (more precisely, satisfy the condition I = /to/2, see (15)). In the strong 
field limit we have 

e£ * 2 eE t 2 

I ~ v t + ~ -, (24) 

ml ml 

\ D ~ h/(k + eE t) ~ h/(eE t) , (25) 
and equating I with /Cd/2 we obtain the "meeting time" t\ and the "meeting coordinate" l±: 

( mh\ 1/3 h h fe 2 E 2 \ 1/3 

1 ~ \e 2 E 2 ) ' 1 " 2eE h ~ 2eE \ mh J ' [ } 

For the energy obtained by the particle from the field in the course of its production we find 

5e D ~ eEoh = ( . (27) 
\ 8m / 

The parameter that governs the modification of the total probability of the process in the 
strong field limit is therefore 



5e D eE^lx ( e 2 Ef,H 



2 P 2^2\ V3 

= X 2/3 ■ (28 



Here the tilde over Ed is to distinguish this quantity in the strong field limit (x> 1) from 
the weak- field value given in eq. (17). 



2.3 Estimates of total probabilities 

We have found that the parameters that describe the effects of external electromagnetic 
fields on the total probabilities of quantum processes are x i n t ne weak field limit and x 2 ^ 3 
in the strong field case. Let us now try to quantify these effects. If the field-free rate of a 
process is Wo(sq), the field-induced increase of the phase-space volume of the process would 
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imply that the rate of the process in the presence of the field is given by W ~ Wo(£o + 6ed)- 
Let us first consider the weak field limit x <C 1, i.e. 5ed <C Eq. In this case we have 

W ~ Wo(e + 5e D ) = W (e ) + W^e )5e D + ^W£(e )(5e D ) 2 + . . . (29) 

The first derivative term, which is linear in the field strength, vanishes upon the averaging 
over the phase of the field corresponding to the particle production time (or over the angle 
between E and the particle momentum k); the same is also true for all odd-order derivative 
terms in (29). Assuming that Wq has a power-law dependence on Eq and taking into account 
that Sed — X £ o> we therefore obtain 

W ~ W {e )[1 + \w^e )x 2 + \w { , IV \e,) X A + ...]■ (30) 

As an example, consider allowed nuclear /3-decay with non-relativistic energy release, such 
as tritium /3-decay. In this case Wq{eq) oc e (see eq. (42) below), and eq. (30) yields 

W * ^o(. )[l + |x 2 + ^X 4 + ...]. (31) 

This expression does not depend on the frequency of the external field, i.e. it actually 
corresponds to the limit u — » 0. How should the dependence on u enter into the expression 
for Wl The parameter governing this dependence is the ratio of the characteristic time 
scale of the process t x and the field period T, i.e. it is ~ ut x . Due to the time reflection 
invariance of QED, the probability of the process can only depend on the even powers of 
this parameter. Thus, the coefficients of \ 2n m t ne expression for W should actually be 
power series in (ut x ) 2 . The coefficients in eq. (31) are just the leading (zero order) terms in 
these expansions. Note that in the considered case x< 1 we have t x ~ to = ^/^Eq, so that 
(ut x ) 2 = (hw/2E ) 2 . 

Consider now the strong field limit )( > 1, in which 5e~d — X 2 ^ £ o- 111 this case to 
leading order one can neglect e in the expression e ~ e + 5ed and expand W (e) in powers 
of e /5e d = x~ 2 ^- This gives 

W ~ W (e + 6e D ) ~ W (e oX 2/3 )[1 + 0( X ~ 2/3 )] . (32) 

For allowed nuclear /3-decay with emission of non-relativistic charged leptons, in which 
Wq(eo) oc e 7 J 2 , we find 

W ~ W (e ) ■ X 7/3 ■ [1 + 0( X ~ 2/3 )} ■ (33) 

Just like the weak- field expression (31), this result actually corresponds to the limit u — > 0. 
The w-dependence of the decay rate W should be given by the power series in the parameter 
(ut x ) 2 , where the characteristic time t x is now given by the "meeting time" ti defined in 
eq. (26), This yields 
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This means that the u dependence of W should only enter starting the term of order x _4//3 
in the expansion (33), whereas the leading term and the term ~ x~ 2 ^ 3 in the square brackets 
should be independent of the frequency of the external field. 

Note that the leading term in (33) (~ Wq{eq)x'^) could actually have been guessed. 
Indeed, since Wq(eq) oc e 7 J 2 and x 7 ^ oc e 7 ^ 2 , the leading term in (33) is independent of 
Eq. This is exactly as it must be in the strong field limit - in this case the field- induced 
correction to the charged particle energy is large compared to the field-free energy, and the 
probability of the process should be approximately independent of the latter. 

It is interesting to note that the fact that the leading-order term in W is independent 
of £o means that is is also independent of its sign. Therefore, eq. (33) holds true even in the 
case Eq < 0, which corresponds to the situation when in the absence of the external field 
the system is stable (or the reaction does not go because the available energy is below the 
threshold). For example, a sufficiently strong field can cause /3-decay of an otherwise stable 
nucleus, such as 3 He. The decay rate for this case is also described by eq. (33) provided 
that by Wo and x one understands WoG^ol) an d x(| £ o|)- In other words, a very strong field 
"pulls" the electron (or positron) out of the nucleus irrespective of whether the nucleus was 
/3-active or stable in the absence of the field. 

2.4 Weak-field limit for systems stable in the absence of the ex- 
ternal field 

Consider now in more detail the case of negative energy release (e < 0), e.g. field- induced 
/3-decay of an otherwise stable nucleus. One can model such a situation by a particle 
bound in the potential well (fig. 2). In the absence of the external field the particle is in a 
stationary bound state. The external dipole electric field adds the potential U = —eE(t)x 
to the potential well, thus transforming it into a potential with a barrier. The state of the 
particle then becomes quasi-stationary; the particle can escape from the well either through 
the tunneling effect or due to a multi-photon "ionization" . 



Figure 2: Left: a particle bound in a potential well; right: a particle in a well with a 
potential barrier. 

Let us start with a qualitative analysis, similar to the one performed in the case Eq > 0. 
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The /3-particle can be emitted (along with the neutrino) virtually, and has to be re-absorbed 
after the time to ~ Vkol- O n the other hand, the time it takes for the electron (or positron) 
to obtain the missing energy |eo| from the external field is t 2 = \/2m\eo\/eE. If this time is 
smaller than or of the order of the "time of virtual existence" of the /3 particle to, i.e. 



t 2 = ^2m\e \/eE < h/\e \, (35) 

the induced decay occurs with a large probability (W/Wo(|£o|) ^ !)• 111 the opposite case 
t 2 3> to the decay probability is strongly suppressed. 

The condition in eq. (35) is equivalent to x(l £ o|) ~ 1 ( see (17)), i.e. corresponds to the 
strong field limit which was considered in the previous subsection. Let us now concentrate 
on the weak field limit x(|eo|) <C 1. in this case the nature of the suppression of the rate 
of the induced /3-decay will depend crucially on the value of ut 2 . For ut 2 <C 1 the external 
field is nearly static, and the field-induced /3-decay proceeds as a tunneling effect. From 
fig. 2 we find that the width of the barrier l 2 satisfies eEol 2 = \eq\, i.e. 

h = \e \/eE . (36) 

The parameter t 2 defined in (35) is then just the tunneling time. The condition x C 1 
means that the WKB approximation can be used, and for the tunneling probability we find 



W ~ exp 



2 

h 



\j2m(\e Q \ - eE x) dx = exp ^ - — J . (37) 



Note that this result is non-perturbative in the field strength Eq. Thus, in the low- frequency 
limit ut 2 1 the probability of the field-induced /3-decay is exponentially suppressed. This 
result applies quite generally to all processes that do not occur in the absence of the external 
field. The pre-exponential factor in the expression for W cannot be found from simple 
considerations like those presented above. The field frequency dependence of W should, as 
usual, be given by the parameter (u>t x ) 2 , where t x now is the tunneling time t 2 defined in 
eq. (35), which gives 

(cot x ) 2 = C 2 - (38) 

Consider now the weak field limit in the case ut 2 1. The potential barrier then 
oscillates fast on the time scale of the tunneling time, and the field-induced /3-decay proceeds 
via the absorption of many photons by the /3-particle, so that this particle goes over the 
barrier instead of tunneling through it. The minimum number of the absorbed photons is 
therefore ri\ = \\e \/ftw~], i.e is the smallest integer that is > \e \/ftw. The probability of 
the process should then depend on the field strength as 

W oc E 2r|£ol/M , (39) 

i.e. it should exhibit a power-law rather than exponential suppression. 
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3 Allowed /3 decay in the field of a strong wave 



We shall now go from simple estimates to direct calculations of the total probabilities 
of quantum processes in a field of a strong electromagnetic wave, taking allowed nuclear 
/3-decay as an example. 

Let us first recall the decay rate calculation in the field-free case, According to the 
Fermi's Golden rule, the decay probability per unit time is given by 

W (e a) = | \M„r / - M lt ? - E e - K) , (40) 

where p and k are the momenta of the emitted neutrino and electron, E e and E e are 
their full energies, and M, and Mf are the masses of the parent and daughter nuclei. The 
energy release in the process is Eq = (Mi — Mf — m)c 2 , where m is the electron mass (we 
neglect the recoil energy of the final-state nucleus). The matrix element of the process 
Mfi = (Gp/\/2) cos6cM.fi, where Gf is the Fermi constant, 6c is the Cabibbo angle, and 
Aifi is the nuclear matrix element (Mfi = 0(1) for allowed /3-transitions). The 5-function 
in (40) allows one to remove one of the energy integrations. Let us recall how it appears in 
eq. (40). The time dependent amplitude of the process is A(t) oc e^ E f~ Ei ^. The squared 
modulus of the integral of the amplitude over the production time therefore yields 

/oo 2 /*oo roo 

dtA(t) oc / dh dt 2 e^- E f ){tl - t2) = T-2tvH5(E 1 - E f ), (41) 
-oo J —oo J — oo 

where we have switched from the integration variables t\ and t 2 to (t\ +^2)/2 and r = t± — 1 2 . 
The integral over (t\ + t 2 )/2 is trivial and yields the normalization time T that has to be 
sent to infinity at the end of the calculations. The decay rate (i.e. the probability per unit 
time) is obtained by dividing the total probability by T and therefore is T-independent. 
The integral over the difference of t\ and t 2 gives the 5-function 5(Ei — Ef), which helps 
perform the phase space integration in eq. (40). 

In the non-relativistic case (eo mc 2 ) the direct calculation of (40) yields 

W^o) = ^(2m) 3 / 2 (4/10 57 r 3 )4 /2 , (42) 

i.e. Wq oc e 7 J 2 , as was pointed out above. 

How can one calculate the decay rate in the presence of an external electromagnetic 
wave? Usually, the calculation is done by replacing the standard field-free wave function of 
the emitted electron with the exact solution of the Dirac equation in the field of a monochro- 
matic electromagnetic wave - the so called Volkov solution [9]. In the non-relativistic case 
and assuming that the dipole approximation is valid, one can use instead a much simpler 
wave function [10], which in the Coulomb gauge (2) it can be written as 

^ k (r,t) = exp|^kr-^| ^ [k-(e/ C )A(t')] 2 rft , }. (43) 
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This is an exact solution of the Schrodinger equation for an electron in a uniform electric field 
of arbitrary strength and arbitrary time dependence. Since in a non-stationary external field 
the ^-dependence of the electron wave function is \l/fc(r, t) oc e~ 1 ^ ^ e~^ Et , the calculation 
of the amplitude does not lead to an energy- conserving 5-function. This simply reflects the 
fact that the energy of a system in a time-dependent external field is not conserved. For a 
periodic electric field E(t) with a period T = 2tt/co the time-dependent factor e" % ^^ in the 
electron wave function can be expanded in a Fourier series: 

oo 

e -</(t) = a n e- inult . (44) 

n=— oo 

Substitution of this into the expression for the transition amplitude would lead to a repre- 
sentation of the amplitude as a sum of partial amplitudes corresponding to the absorption 
from the field or stimulated emission into it of all possible numbers of photons. 

Consider a circularly polarized wave with the electric field strength 

E(t) = {E smut, - E cosut, 0} (45) 
By making use of the above Fourier expansion technique one can present the decay rate as 

- - >^S±j @) *• + *-*-.-« • (46) 

where k± = ^k 2 . + fc| is the component of the electron momentum in the plane of the field, 
e = k 2 /2m is the kinetic energy of the electron, Ek = e^g^mw 2 is the mean kinetic energy 
of the oscillatory motion of the electron in the field of the wave, and n = [(e — Ek)/Hu\, 
i.e. the integer part of (eq — Ex)/hu. The appearance of the lower limit —no in the sum in 
(46) is related to the fact that the electron cannot emit into the field more energy than it 
has. In the limit Eq — > only the term with n = in the sum survives and, since Jo(0) = 1, 
the field-free probability (40) is recovered. 

The expression in eq. (46) was obtained by making use of the standard technique of 
Fourier expansion of the electron wave function. Now, let us look at this formula more 
closely. It contains an integral over the electron and neutrino 3-momenta of the sum of 
squared Bessel functions, whose arguments depend on the electron momentum and on the 
angle between this momentum and the direction of the external field strength. Moreover, 
because of the presence of the (^-functions, the argument of the Bessel functions J n also 
implicitly depends on the index n (through the dependence of k on n). The summation 
extends from a finite value of n to infinity. Obviously, it would be very difficult to calculate 
the probability of the process using this expression! 

It is actually easy to understand why this calculational nightmare arises. The reason 
for this is a bad calculational approach. Indeed, each term in the sum in (46) gives a 
partial probability of /3-decay with absorption from the field (for n > 0) or emission into it 
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(for n < 0) of \n\ photons. The energy conservation for each such processes is ensured by 
the corresponding ^-function in the integrand. In fact, this means that the calculation is 
done in the energy representation. However, for a system in a non-stationary external field 
energy is not a good quantum number, therefore using the energy representation does not 
give us any advantage. Moreover, in order to calculate the total probability of the process, 
it is actually not necessary to calculate all partial probabilities and then sum them; the 
partial probabilities actually contain much more information than is needed. Equation (46) 
therefore uses excessive information. 

The above arguments actually give us a hint of how an adequate calculation should be 
done. One should not make use of the Fourier expansion and calculate the partial probabil- 
ities. There will be no energy- conserving ^-functions in this case, but the integrations over 
the electron and neutrino momenta can be performed without using 5-functions. The idea 
of a simple calculation presented below is the following: 

1. In the calculation of the squared modulus of the integral of the amplitude of the 
process over the production time go from the variables ti and t 2 to it\ + t 2 )/2 and 
t — ti — t 2 ? & s i n the calculation that led to eq. (41). 

2. Perform first the integration over (ti + t2)/2, then over the electron and neutrino 
3-momenta, and only at the very end do the integration over the difference of the 
times r. 

Unlike in the field- free case, the integral over (ti + 1 2 )/2 is not trivial now, but it still leads 
to the proportionality of the total probability of the process to the normalization time T. 
For the decay rate (i.e. probability per unit time) we find 

ft 2 1 N J {2nhy J {2nhf 

T / e Eq kj_ . ut\ r i ( k 2 e 2 El _, \ 1 , , inS 



Note that at this point it is still easy to go to the standard calculational technique. Indeed, 
by using the relation [11] 

oo 

J (2xsmy)= £ J 2 (x)e 2in y (48) 

n=— oo 

one can reduce eq. (47) to (46). 

We now perform in (47) the momentum integrations, leaving the integration over r = 
t\ — £2 to the end. The integration over the neutrino 3-momentum is trivial, 2 whereas the 



2 In this calculation the neutrino mass is neglected. We discuss the case m„ 7^ in sec. 3.4. 
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integral over the electron momentum can be readily done in cylindrical coordinates. The 
result is [21 



W 



m 



3/2 \M fi \ 2 (h 



1 7/2 



2 5 ir A H 7 c 3 

Here the integration variable is x 
parameters 

6-. 



dx 



x (x + a)*/* exp 1 " 1 r + - x )\\- (49) 

u(ti — £2)/2, and we have introduced dimensionless 



2£o 



7 



e 2 El 
mhu 3 



(50) 



Note that in terms of 7 and 5 the parameters £ and x introduced in sec. 2 are expressed as 



^ s 



x 2 = i/s 3 



(51) 



Equation (49) is the result we were looking for. Instead of the frightening expression (46) 
we have now obtained a relatively simple single integral. It can be calculated numerically, 
but all the interesting limiting cases can actually be readily studied analytically. 

Let us first note that our approximation (1) implies <5 3> 1; therefore, the integral in 
(49) gets its main contribution from the region |x| < <5~ x 1 and from the vicinities of the 
stationary phase points which lie outside this region. The stationary phase contributions 
are strongly suppressed fast-oscillating functions of the field strength Eq which vanish upon 
averaging over small fluctuations of this quantity; we will discuss their contribution later 
on. For small |x| one can expand sin 2 x in the exponent in eq. (49) in powers of x. Keeping 
the first two terms in this expansion, we find 



W 



m 



3/2 \M fl \ 2 (hi 



2 5 7T 4 ft 7 C 3 

The integral here is of Airy type. 



.7/2 



dx 



[x 



iO) 9 / 2 



cxp 



{-,-[ < 5.x- 7 y]}. (52) 



3.1 Weak-field limit (case £q > 0) 

Consider first the limit 

7 « 5 3 (x « 1) • (53) 

For brevity, we shall call it 'the weak-field limit' even though the field- dependent parameter 
7 may actually be much greater than unity. For the values of x satisfying \x\S < 1 one then 
has 7|x| 3 1. Therefore in this regime one can expand 



cxp 



x 



S-x — 7 — 

' 3 



iSx 



X 



1 + z 7 — + . . . 



Substituting this into (52), we find 

oc 

W = W (e ) ■ 



r(9/2) 



,2m 



m=0 



m!3 m r(9/2 - 3m) 



X 



(54) 



(55) 
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The first few terms in this expansion give 



W 



W (e ) 



1 + 



35 



-X 



35_ 4 
128 X 



(56) 



Comparing this with our estimates made in sec. 2 (see eq. (31)), we find a very good 
agreement. Not only the facts that the decay rate W depends on the characteristics of the 
external field through the parameter x an d that for small x it is a power series in x 2 were 
predicted correctly by our simple estimates - even the numerical values of the coefficients 
of the x 2 and and x 4 terms were found correctly. This gives an additional a posteriori 
justification to our choice of the numerical coefficient in the expression for the characteristic 
length l x in eq. (15). 

The expression in eq. (56) actually correspond to the stationary field limit uj — > 0. To 
obtain the dependence of W on u one has to retain next terms in the expansion of sin 2 x in 
the exponent in eq. (49). For the lowest order terms we obtain 



35 

W = W (e ) I I + —x 



i - if^V 

30 \2e J 



+ o( x 4 



(57) 



i.e. in the weak field limit the (^-dependence enters into the expression for W as an expansion 
in powers of (hiu/2eo) 2 , as it was anticipated in sec. 2. 

Let us now turn to the contributions to W from the stationary phase points. From (52) 
one founds two stationary phase points: 

±(5/7) 1/2 



= ±£ 1 . 

Their contribution to the decay rate is, to leading order, 

105 

SW = W (e )— X 4 cos(2/3 X ) 
16 



(5? 



(59) 



Thus, the stationary phase contributions are suppressed at least as x 4 and are fast-oscillating 
functions of the field strength Eq with zero average. 



3.2 Strong field limit (x > 1) 



In the strong field limit X 3> 1 (which corresponds to 7 ^> 5 3 ) one can, instead of (54), use 
the expansion 



exp 



5-x — 7 



e' 7 T 



i5x 



5 1 
2 



Substituting this into eq. (52) yields 

35 3 5 / 6 ^ (_ 3 )m x -2m/3 
m=0 



w 



ml 



2nm n 



in 



(60) 



(61) 
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This result holds for /3-active nuclei (eq > 0). For nuclei that are stable in the absence of 
the field (e < 0) the decay rate can be obtained from eq. (61) by replacing (— 3) m — > 3 m , 
W (e ) -> W (\e \) and x(e ) x(M)- 

To obtain the dependence of the decay rate on the field frequency u one has, as usual, 
to retain the next terms in the expansion of sin 2 x in powers of x in eq. (49). This yields 



^ = ^o(kol)x 7/3 



7T 



+ 



35^ 
^8" 



-X 



-4/3 



5-3 5 / 6 -r(5/6) 7-3 1 / 6 -r(l/6) _ 2/3 

16v^F X 



+ 0(x 



1 + 

15 \2soJ 



(62) 



Here the upper and lower signs correspond to Eq > and Eq < 0, respectively. 

Let us compare this result with the estimates made in sec. 2. It was predicted there 
that in the strong field limit W ~ W (e )x 7 ^ 3 times a power series in x~ 2 ^ 3 ( see (33)). 
Equation (62) confirms this result, except that the leading term has an extra constant 
factor C\. However, for the numerical value of this factor we have 



5-3 5 / 6 -r(5/6) 



0.9943 



7T 



(63) 



which is very close to 1. Thus, the predictions based on simple estimates are confirmed 
once again. The first two terms in the curly brackets in eq. (62) are w-independent, and the 
dependence of W on u> comes through the positive powers of the quantity x~ 4//3 (^ w /2£o) 2 ) 
again in full agreement with the results of our qualitative analysis in sec. 2. As was dis- 
cussed there in detail, the strangely-looking dependence of W on fractional powers of x 
is a consequence of the fact that in strong fields the electron's de Brogle wavelength is 
field-dependent. 



3.3 The case of weak fields and £q < 

Consider now the case of relatively weak external fields and nuclei stable in the absence of 
the field (eq < 0). As was discussed in sec. 2, there are two essentially different regimes of 
the induced /3-decay in this case, the tunneling regime and the multi-photon one, depending 
on whether the field frequency u is small or large compared to the inverse tunneling time 



t^ 1 = eEo/ y/2m\£o\. Both cases can be obtained from the same master equation (49). In 
the low-frequency limit ut2 = <C 1 in the zeroth order in this parameter one can actually 
make use of a simpler formula (52). The calculation of the integral in (52) in the steepest 
descent approximation yields 

k — ft — 
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To take into account the u dependence of the decay rate, one has to retain the higher-order 
terms in the expansion of sin 2 x in the exponent in eq. (49). This gives, to leading order, 

W(x) = W (\e \)—(l + — 2 ) x 4 exp 

In the multi-photon regime (wia = £ _1 ^* 1), one can find the probability of the process 
by calculating the integral in (49) in the saddle point approximation. The simpler expression 
(52) cannot be used in this case since a region of complex x with large modulus is important 
in the integral. The saddle point is found from the transcendental equation 

sinh 2 z — (cosh z — z" 1 sinh z) 2 = £ -2 , (66) 

where z = ix. Solving it approximately in the limit ^ < 1, we find 

Txrn ,,105 f frw V f, 2m\e \uY 9/2 ( eE \ ^" 1/2 

Thus, we see that, in accordance with the qualitative analysis of sec. 2, in the weak-field 
limit x ^ 1 the probability of the field-induced (3 decay is exponentially suppressed in the 
low- frequency case and exhibits a power- law suppression in the high-frequency limit. 



— ( 1 

'3x V W 



(65) 



3.4 Neutrino mass effects 

Up to now, in our calculations we have been neglecting the neutrino mass m v . In the 
absence of external fields, the only characteristic of /3-decay with the dimension of energy 
is the energy release Eo, and it is always much greater than m„c 2 . This justifies the neglect 
of m v in calculations of the total /3-decay rates in the absence of external fields. 

In the case of /9-decay in a field of monochromatic wave, one can construct other quan- 
tities with the dimension of energy, which depend on the field strength and frequency and 
which have to be compared to m u c 2 . If any characteristics of the process should turn out 
to depend on the ratios of m u c 2 and these field-dependent parameters, study of /3-decay 
in electromagnetic fields could provide important information on the value of the neutrino 
mass. 3 Using E v = (p 2 + m 2 ) 1 / 2 in eq. (47) and following the same steps as those that led 
to (49), we obtain [2] 

w = -^^^ 3/2 \MM 2 K,c 2 ) 2 (M 3/2 

f°° dx ( m u c 2 \ f r /sin 2 x \"n . . 

x I. ¥TW 2 Y 2l ^r x ) exp i " 1 r + H— - x ) \ I ■ (68) 

3 Note that the effective neutrino mass m v that can be probed in /3-decay is expressed through the masses 
of neutrino mass eigenstates m, and the elements of the leptonic mixing matrix U as = \U e i\ 2 mi. 



21 



This equation can be studied by the same methods as those that were applied to eqs. (49) 
and (52). The analysis [3] shows that the dependence of the /3-decay rate W on the neutrino 
mass has the same nature as its dependence on the field frequency u, i.e. it enters through 
the parameter (m u c 2 t x /h) 2 , where t x is the characteristic time scale of the process. As was 
discussed above, t x is different in different regimes. The parameter (m u c 2 t x /h) 2 turns out to 
be very small in all cases except for e < 0, x ^ 1) when t x is the tunneling time £2 defined 
in (35). In this case the corrections to the probabilities of field-induced /3-decay due to non- 
vanishing neutrino mass can be sizable. However, the decay probabilities themselves are 
extremely small then, making the neutrino mass effects unobservable. Thus, unfortunately 
/3-decay in external electromagnetic fields cannot tell us much about the neutrino mass, at 
least as far as the total probabilities are concerned. Note that interesting effects on the 
spectra of /3-electrons may still be possible. 

4 Forbidden (3 decay in strong fields — lifting the for- 
biddenes 

As was pointed out in sec. 2.2, one possible mechanism of modification of the total probabili- 
ties of quantum processes in electromagnetic fields is modification of their transition matrix 
elements. We shall discuss now a particular example of this - acceleration of forbidden 
nuclear /3-decay due to the interaction of the parent or daughter nuclei with the field. 

When considering allowed /3 decay in external fields we were taking into account only 
the interaction of the produced electron or positron with the field; the interaction of the 
involved nuclei was ignored because they are very heavy. In the case of forbidden /3-decay, 
however, there exists an enhancement mechanism which requires taking the interaction of 
nuclei with the field into account. If the parent or daughter nucleus absorbs from the field 
(or emits into it) one or more photons in the course of /3-decay, the angular momentum of 
the initial or final nuclear state may change, and the forbiddeness may be lifted. 

Let us discuss this mechanism in more detail. Let a nucleus undergo a forbidden /3-decay 
from its ground state \i) to the ground state |/) of the daughter nucleus. Assume that 
the parent nucleus has an excited state |1) whose quantum numbers permit an allowed 
/3-transition |1) — > \f). The parent nucleus can then undergo a virtual transition to the 
state |1) by absorbing one or more photons from the external field (or by emitting them 
into the field), followed by the allowed /3-transition. An analogous situation will obtain if 
the daughter nucleus has an excited state |2) with quantum numbers permitting an allowed 
/9-decay \i) — > |2), followed by the electromagnetic |2) — > |/) transition under the influence 
of the external field. 4 

4 We are assuming here that the excitation energy of the state |2) is larger than the energy release eo 
of the ground state — » ground state /3-transition. Otherwise the allowed /3-decay \i) — > \2) followed by a 
cascade of 7-transitions |2) — > |/) would be possible even in the absence of the external field. 
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Let us now concentrate on the simple case of unique first-forbidden /3-transitions for 
which the selection rules are AJ An = 2~. Higher-order forbidden transitions can be exam- 
ined similarly. 

For unique first-forbidden transitions, the state |1) must be related to the ground state 
\i) of the parent nucleus by the electric dipole (El) transition. The |1) — > \ f) transition 
will then be the allowed Gamow- Teller transition with AJ Aw = 1 + . An analogous situation 
would arise for transitions through the virtual state |2) of the daughter nucleus. It is 
readily seen that any other possibility (e.g., electromagnetic M2 transition and allowed 
Fermi transition with AJ An = + ) would result in a much smaller matrix element. 

Let us estimate the matrix element Mi n d of the induced process corresponding to the 
relaxation of the forbiddeness in external fields. We shall denote the energy differences 
between the states |1), \i) and |2), |/) as 

Aei = £1 - Ei , Ae 2 =e 2 — e f . (69) 

The admixture of the state |1) to the ground state \i) in the external field of frequency 
oj <C Ae\ is characterized by the parameter eEodu/ Ae\, where du is the dipole matrix 
element of the corresponding electromagnetic transition. Similarly, the admixture of the 
state 1 2) to the ground state |/) of the daughter nucleus is characterized by the parameter 
eE d 2 f I ' Ae 2 . To estimate the matrix element M in ^ we replace the dipole matrix elements 
d\i and d 2 f by the nuclear radius R and obtain 

M ind ~ ^ • 1 , (70) 

where for the estimate we have assumed that the matrix elements of allowed /3-transitions 
are of order unity (in reality, for intermediate and heavy nuclei they are somewhat smaller) . 

Note that the "forbiddeness lifting parameter" eE R/ 'Ae 12 does not depend on the field 
frequency oj in the limit oj <C A£ 12 . This is in contrast with the results of refs. [12, 13], 
where it was claimed that the forbiddeness can be relaxed if the parameter z 1 / 2 = eE^Rjoj 
becomes sizeable. As was discussed in sec. 2, the dependence of the total probabilities on 
the parameters that diverge in the limit oj — > is not admissible from the physical point 
of view; it was shown in [3] that the results of refs. [12, 13] were a consequence of some 
unjustified approximations adopted in those papers, which, in particular, led to a breakdown 
of gauge invariance. As a result, the probabilities of forbidden nuclear /3-decay in external 
electromagnetic fields were overestimated in [12, 13] by many orders of magnitude. 

The dependence of the induced matrix element Mi n d on the external field frequency can 
become important for large enough oj, when u ~ Aei >2 \ in this case in eq. (70) Ae 12 in 
the denominator has to be replaced by Ae 12 ± oj. In particular, for oj = Ae 12 a resonant 
enhancement of the process becomes possible; the denominator of M ind in eq. (70) would 
then just contain the width of the corresponding excited state, Ti or T 2 . Unfortunately, 
such a resonant enhancement would require 7-lasers, which do not exist at present. 
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Let us estimate the enhancement factor due to the relaxation of the forbiddeness of 
the /3-transition in low-frequency fields. In the absence of the field, the matrix element of 
unique first-forbidden /3-decay is M ~ koR/h, where, as usual, ko = \/2meQ. From eq. (70) 
we then obtain 

M ind ^ eE h _ 

M ~ y/2F^Ae lt2 ~ X *' 1 ' 

i.e. the parameter that governs the enhancement of the probability due to the lifting of the 
forbiddeness is x* — x(% £ o/ Aei,?)- For the rate of the process in the case of relatively weak 
fields we then expect 

W ~W (£ )[l + a X 2 + bxt + ...}, (72) 

where a and b are numerical constants. 

Let us now sketch the calculation of the decay rate (the details can be found in ref. [3]). 
Although all physical observables are gauge invariant, a judicious choice of the gauge can 
simplify the calculations greatly. For calculating the probability of forbidden /3-decay in 
external electromagnetic fields, the scalar gauge (3) turns out that to be most convenient. 
The reason for this is that the corrections to the nuclear wave functions can then be described 
in the lowest order in the small parameter eE R/ 'Aei2 (see (70)). At the same time, the 
corresponding interaction parameter in the Coulomb gauge, 5 



e A p ni 1 _ eEpdn j 
c M (e„ — Si) huj 



(73) 



is not in general small in low-frequency fields, and therefore it cannot be used as an expansion 
parameter in perturbation-theory calculations. (Here Aq = cEq/uj is the amplitude of the 
vector-potential and n corresponds to an excited nuclear state). Moreover, unlike the scalar- 
gauge parameter eE R/ Ae± t 2, it is not suppressed for large values of e n — Ef, this means that 
in calculating nuclear wave functions in the presence of external fields one cannot constrain 
oneself to the contributions of only lowest-lying excited states and has to sum over the whole 
spectrum of nuclear excitations. Therefore, we choose to calculate in the scalar gauge (3). 

The wave function of a non-relativistic electron in the scalar gauge takes the form [10] 

K(r, t) = exp {^[k - (e/c)A(t)]r - ^ j\k - (e/c) A(t')] 2 df } , (74) 

where the prime refers to the quantities in the scalar gauge. Note that here A(t) does not 
have the meaning of a vector-potential, which is zero in the scalar gauge (3); by definition, 
in eq. (74) A(t) = -c fdfEtf). 

Comparing the electron wave functions (43) and (74), we find that they are related by 

*i(r I t) = e-fe A ^ fc (r,*). (75) 



Here we consider the parameter describing the modification of the wave function of the parent nucleus. 
The daughter nucleus can be considered similarly. 
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This can be readily understood by noting that the gauge transformation from the Coulomb 
gauge (2) to the scalar gauge (3) is A^(t, x) — > x) + d^r](t, x) with the gauge function 
r)(t, x) = A(t)r. The wave functions of particles of charge e should then transform as 

V>(t, r) e^ ie/hc)v{t ' r) ij(t, r) . (76) 

The electron wave functions in the Coulomb gauge and scalar gauge thus have different 
coordinate dependence. For allowed /3-transitions the electron wave function inside the 
nucleus is replaced by a constant, and therefore it does not matter whether the Coulomb- 
gauge or scalar-gauge expression is used. For calculating the rate of forbidden /3-decays the 
wave function (74) should be employed. 

The calculations [3] confirms the expected result (72) for the rate of unique first- 
forbidden (3 decay in the field of a strong electromagnetic wave. Since nuclear excitation 
energies Ae 12 are typically of the same order of magnitude as the energy release e , and 
in all known cases are not smaller than ~ 1 keV, for currently attainable fields the forbid- 
deness lifting term bxl in (72) is extremely small, as is the phase space enhancement term 
ax 2 - We discuss this point in more detail in the next section. Notice, however, that the 
numerical value of the coefficient a in (72) is significantly larger than in the case of allowed 
/3-transitions: for unique first-forbidden /3-decays the calculation [3] yields a = 315/8, to be 
compared with 35/8 in eq. (56). 

5 Summary and discussion 

We have considered, both qualitatively and quantitatively, the influence of the field of a 
strong electromagnetic wave on the characteristics of quantum processes with participation 
of non-relativistic charged particles. Our qualitative analysis has shown that the parameter 
that determines the influence of the external field on the differential characteristics of the 
processes (such as energy spectra and angular distributions of the final-state particles) is 
£ = eE / '(i/2me uj). The processes of modification of these differential characteristics in 
external fields are essentially of classical nature. 

At the same time, the parameter that governs the modification of the total probabilities 
of the processes (such as decay rates and total cross sections) is x = gE q H/ (i/2me 2e )- 
It describes the energy 5ed obtained by the charged particle from the field (or given to 
the field) over the distances of order of the particle's de Broglie wavelength. In relatively 
weak fields, when ^ < 1, this energy is given by 5eo — X £ o, whereas in the strong field 
limit (x ^> 1) one has 5e~d — X 2 ^ 3£ o- The fractional power of x enters in the latter case 
because of the field dependence of the charged particle's de Broglie wavelength. The process 
of energy exchange between the produced particle and the external field in the process of 
the formation of the particle (and therefore the modification of the total probability of the 
process) is inherently quantum in its nature. 
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We also estimated the total probabilities of quantum processes in external electromag- 
netic fields, considering nuclear /3-decay as an example. Our estimates were made in the 
limits x < 1 and % ^ 1- I* 1 the latter case the decay rate for a /3-active nucleus (eq > 0) 
and the rate of the field-induced /3-decay of the daughter nucleus which is stable in the 
absence of the external field (eq < 0) are practically the same. Examples of such transitions 
are 3 H ->■ 3 He + e" + u e (e > 0) and 3 He ->• 3 H + e + + u e (e < 0). In the weak- field 
limit the induced decay in the case Eq < is predicted to be exponentially suppressed in 
the low-frequency (tunneling) limit ut 2 <C 1, where t 2 is the tunneling time defined in (35), 
and power-suppressed in the multi-photon regime ut 2 3> 1. 

The dependence of the total probabilities of quantum processes on the field frequency u 
comes through the expansion in powers of the parameter (ut x ) 2 , where t x is the characteristic 
formation time of the produced charged particle. At the same time, as was mentioned above, 
the energy obtained by the particle in the process of its formation is given by the work of 
the field on the particle over the distance of order of the particle's formation length l x : 
5e ~ eE l x . The quantities l x and t x are different in different regimes: 

• Weak fields : 

E > : l x ~ J= , t x ~h/2E ut x = huj/2E = 1/5, 

2^2m£ 



kol \/^ m \^o\ i i — r 

e < : l x ~—=-, — f =^ w ^ = V 2me ow/e-Eo = ■ 

eE eE 

• Strong fields : 

/ h 2 \ 1/3 / rah \ 1/3 

We have also discussed a simple method of calculating the total probabilities of quantum 
processes with participation of non-relativistic charged particles in the field of an electro- 
magnetic wave, considering nuclear /3-decay as an example. The method does not rely on a 
summation of partial probabilities with absorption from the field or emission into it of all 
possible numbers of photons; instead, the total probabilities are calculated directly. The re- 
sults of the direct calculations in the cases of allowed and forbidden /9-decay fully confirmed 
the estimates made in sec. 2, often even including the values of numerical coefficients. 

How large can actually the effects of strong electromagnetic fields on /3-decay rates be? 
The parameter x that determines the modification of the rates can be written as 

eE h Eq ( mc 2 ^ 

x = 7' 

where 



X V2rnT 2E E c \ 2e q ' 



m 2 c 3 



E c = —— = 1.323 x 10 ib V/cm (71 
eh 
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is the QED critical field strength (the so-called Schwinger field). The most powerful present- 
day lasers can reach the field intensities up to I ~ 10 22 W/cm 2 , and in the near future 
probably the intensities I ~ 10 24 W/cm 2 will become available. From the formula 

£ (V/cm) ^ 20 v//(W/cm 2 ) (79) 

we then find that this corresponds to the field strengths E /E c ~ 10 -4 - 10~ 3 . For tritium 
/3-decay from eq. (77) we then obtain 

X ~ 5 x 10~ 3 - 5 x KT 2 . (80) 

Taking into account that for small x the correction to the decay rate is of order x 2 , we see 
that the field effect on tritium /3-decay is too small to be observable. 

Since x 2 scales as Eq 3 , one can expect significantly stronger effects for /3-decays with 
smaller energy release. In particular, for unique first forbidden /3-decay 187 Re(| + ) — > 

187 0s(tj ) one has Eo ~ 2.64 keV, and the parameter \ can formally take values x ~ 1 
even for present-day lasers. Does that mean that we can already observe strong effects of 
the laser fields on nuclear /3-decay? 

Unfortunately, in reality the situation is not that promising at the moment. First, for 
nuclei with small energy release (and especially for forbidden /3-decays) the decay rates are 
extremely small, and even sizeable corrections to them are not easily observable: the life- 
time of 18 Re, for example, is ~ 5 x 10 10 yr. The situation is complicated by the fact that 
the powerful lasers have very low pulse duration (in the femtosecond range) and repetition 
frequency which is typically ~ 10~ 3 Hz. Finally, except for oj 3> u at where u at are char- 
acteristic atomic frequencies, the atomic electrons would screen the external fields, greatly 
reducing their strengths at the nucleus, and to observe the field effect one would first have 
to ionize the atom of the /3-active element. 

In contrast to this, laser field effects on energy spectra and angular distributions of 
electrons and positrons emitted in nuclear /3-decay can be quite significant. In addition, 
strong electromagnetic fields can still influence sizeably the total probabilities of atomic and 
molecular processes with emission of low-energy electrons. An example of such a process is 
a photo-ionization in two fields, one of which is a weak field with the energy of the quantum 
Ml slightly above (or slightly below) the ionization potential / and the other - an intense 
field with the frequency satisfying ftw I. 
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Appendix A: Validity of approximations 



We consider here the validity conditions for the approximations adopted in our calculations 
of the total probabilities of quantum processes. 

The dipole approximation implies that the wavelength of the external electromagnetic 
field is assumed to be large compared to the characteristic length parameters l x involved in 
the problem: A = 2irc/u ^> l x . Let us consider this condition in various regimes. 



(a) x ^ 1) e o > 0: l x — A~d/2 = H/ (2^2meo). The dipole approximation requires 

fc / 2\ 1/2 



2£o «H^J >>1 - (A1) 

which is certainly satisfied for hid < Eq. 

(b) x ^ 1) £ o < 0: l x = \eo\/gEo. The dipole approximation is valid provided that 

<< ^(^]'\ (A2) 



2|e | \2|e | 
(c) x > 1) > or < 0: l x = [k 2 / (8meE )] 1/3 . The validity condition is 

ftw , /o / mc 2 N 



which is satisfied with a large margin. 

Next, we discuss the non-relativistic approximation. In the weak- field limit (% <C 1) the 
standard condition is 

e < mc 2 . (A4) 

However, in the strong field limit x ^ 1 one also nas to make sure that the energy obtained 
from the field by the charged particle during its formation does not take it out of the non- 
relativistic domain: eE l x <C mc 2 . Taking into account that in the strong- field regime 
lx = [^ 2 /(8mei?o)] 1//3 , we arrive at the condition 

E < E c = 1.323 x 10 16 V/cm. (A5) 

In our calculations of the rates of nuclear /3-decay in external fields we were neglecting 
the final-state interaction of the produced electron or positron with the Coulomb field of 
the nucleus. This is permissible when the Coulomb energy Ze 2 /l x = aZhc/l x is small com- 
pared to the characteristic energy e x of the process. Consider now the Coulomb parameter 
aZhc/l x e x in various regimes. 

(a) x <^ 1 ; e > 0: l x ~ h/(2 v / 2me ), e x ~ e . In this case 

aZhc/e x l x ~ AaZc/vo = AaZc(m/2e ) l/2 , (A6) 
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which is (up to the factor 4) just the standard Coulomb parameter. As an example, for 
/3-decay of 3 H we have Z = 2, (m^o) 1 / 2 ~ 3.7, aZhc/e x l x ~ 0.2. 

(b) x<l,£o<0: ^ = |e |/eJ5, e x = |e |, 



aZhc/e x l x ~ 4aZc(m/2|e |) 1/2 x. 

(c) x»l,e >0 or <0: = [/i7(8me£ )] 1/3 , £ x = X 2/3 ko| 

aZhc/e x l x ~ 4aZ(m/2|e |) 1/2 X" 1/3 - 



(A7) 



(8) 



Appendix B: Some useful formulas 

Here we collect some formulas which have been used in calculation done in sees. 3 and 4. 
An integral representation for Jq(z): 



Uz) 



Gradshteyn & Ryzhik [11], 6.631(6): 



x u+1 e ±tax J v (0x) dx 



[a>0, - 1 < Rez/ < 1/2, /3 > 0]. 
Gradshteyn & Ryzhik [11], 3.382.7: 



7T./o 



'd0. 



/3 l 



2a I/+1 



exp 



. v + 1 /3 2 

±Z 7T 

2 4a 



(/3 - %X) V 



dx 



27rp t/ ~ 1 e- 

r(i/) 
0, 



p > 
p < 



(Bl) 



(B2) 



(B3) 



[Rez/ > 0, Re/3 > 0]. 

For an electromagnetic wave with the electric field strength E(t) = {Eq sin cut, — Eq cos tut, } 
the electron wave function in the Coulomb gauge (43) is 



* k (r, t) = exp { ^kr - ^ 



k 2 e 2 Al\ eA . , ' 

t ( /c^ sm uot — k v cos cut) 

2mc 2 / 



2 

2m 



mcuj 



(B4) 



where A = Eqc/uj. 
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